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Oh, Abstract Possible bulk compositions of the super-Earth exoplanets, CoRoT-7b, Kepler-9d, 



o 



o 



and Kepler-lOb are investigated by applying a commonly used silicate and a non-standard 



C/3 ' carbon model. Their internal structures are deduced using the suitable equation of state of 

the materials. The degeneracy problems of their compositions can be partly overcome, based 

, on the fact that all three planets are extremely close to their host stars. By analyzing the 
> ■ 

■ numerical results, we conclude: 1) The iron core of CoRoT-7b is not more than 27% of its 

ly-^ ' total mass within 1 a mass-radius error bars, so an Earth-like composition is less likely, but 



its carbon rich model can be compatible with an Earth-like core/mantle mass fraction; 2) 
, Kepler- 1 Ob is more likely with a Mercury-like composition, its old age implies that its high 

iron content may be a result of strong solar wind or giant impact; 3) the transiting-only super- 
Earth Kepler-9d is also discussed. Combining its possible composition with the formation 
' theory, we can place some constraints on its mass and bulk composition. 
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1 INTRODUCTION 

Rocky planets (like the Earth in our solar system) located in the habitable zone of their stars are presently 
the best candidates for harbouring extra-terrestrial life (Leger et al. 201 1). Therefore, the search for rocky 
exoplanets plays an important role in the present exoplanet detections. Recent observations of exoplanets 



have revealed 758 exoplanets (http://exoplanet.eu/ updated to 03 November 2011), mainly through stel- 
lar radial velocity measurements or photometric detection of planets transiting their host stars. Among 
them, more than 30 exoplanets are with minimum mass < 10 Earth masses (Mq). They are called 'super- 
Earths' and a faction of them may possibly be rocky planets.Recently, the Kepler mission revealed 1235 
planetary candidates with transit-like signatures detected during its first 4-month operation (Borucki et 
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al. 1201 II ). Among them, 68 candidates of approximately Earth-sizes (Rp < 1.25 i?©), 288 super-Earth 
sizes (1.25 i?® < Rp < 2 i?©). Beyond all doubt, we are entering the era of small exoplanets. Some of 
them are expected to be terrestrial in nature. 

The understanding of super-Earths is of great interest to us due to their possible rocky nature. What is 
the bulk composition of these new planets? How were they formed? These are all interesting yet puzzling 
questions. However, the information presently available based on observations is still Umited for under- 
standing their composition. Through transit and radial velocity detections, planetary radii and masses are 
the only windows presently to investigate their composition for most exoplanets. Among the detected super- 
Earths, their masses and the radii show quite different signatures when compared to terrestrial planet of our 
solar system. Their low average density tells us that they must primarily contain some light element which 
is lighter than water They may have a significant gas envelope like Jupiter or Neptune in our solar system. 
However, super-Earths lacking of gas envelops, or 'rocky super-Earths', may be also present in some cases, 
as we will discuss below. 

CoRoT-7b is the first super-Earth with a measured radius, R ~ 1.68 ± 0.09 i?© (Leger et al. l2009l l. Its 
obtained mass combined with radial velocity measurements is M = 4.8 ± 0.8 M© (Queloz et al. |2009l l. 
After 2009, the initial results have been revised by other authors. For example, after revision of CoRoT-7 
stellar parameter by Bruntt et al. (1201 OK it yields M ~ 5.2±0.8 A/©. Hatzes et al. ( 1201 11 1 performed another 
analysis of the data and got a larger mass, M = 7.0 ± 0.5 M©. The works done by Boisse et al. (1201 11 1 
and Ferraz-Mello et al. (1201 11 1 got different results. In this paper, we do not consider those masses and 
use the initial data of CoRoT-7b only. Recently, two transiting super-Earths, Kepler-9d (Torres et al. 1201 1) 
and Kepler-lOb (Batalha et al. [20TT]) were discovered. Kepler-9d was first discovered in 2010 as a planet 
candidate named KOI-377.03 (Holman et al. l2010l) . It was validated as a super-Earth and the reported radius 
is i? = 1.64^Q i?© (Torres et al. 1201 It . Although current spectroscopic observations are yet insufficient 
to establish its mass, a possible mass range can be found by Holman et al. (120101) . The upper mass limit of 
M — 7 M© corresponds to the maximum mantle-stripping limit (Marcus et al. 120 101 ) for a maximally iron- 
rich super-Earth. The lower mass limit is ^3.5 il/© as a volatile-poor rocky planet with a Ganymede-like 
Fe/Si ratio. So we consider all possible compositions of Kepler-9d using this mass range in our paper In this 
paper, we use Kepler-9d as an example to show what can be inferred for 'transiting-only' hot super-Earths. 
Kepler- 10b was discovered in early 2011, the reported mass and radius are Mp = 4.56+i;29 and 
Rp = 1.416!to'o36 -^ffi (Batalha et al. 1201 11 1, respectively. The main parameters of these three exoplanets 
are summarized in Table 1 . 

The bulk composition of an exoplanet can not be uniquely determined by the measured mass and ra- 
dius. When more than two chemical constituents are taken into consideration, a measured radius and mass 
correspond to more than one possible bulk composition, in other words, it is a degeneracy problem. The 
degeneracy problems of bulk compositions of low-mass exoplanets can be partly overcome if they are very 
close to their host stars because some chemical constituents may be ruled out a priori. Valencia et al. (2010) 
discussed the evolution of close-in low-mass planets like CoRoT-7b. Due to the intense stellar irradiation 
and its small size, it is unlikely to possess an envelope of hydrogen and helium of more than 1/10000 of 
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Table 1 Planet Parameters for Kepler-9d (Torres et al. EoTll . Kepler- 10b (Batalha et al. lSoTTT l. 
and CoRoT-7b (Leger et al. |20091 l. The Ages of Their Host Stars are Also Listed Here. 



Planet Name 



Kepler- 10b 



Kepler-9d 



CoRoT-7b 



Mass (M^) 
Radius (J?®) 
Orbital Period (days) 
Equilibrium Temperature (K) 
Orbital Semimajor Axis (AU) 
Host Star Age (Ga) 



+0.000004 



0.83 7495 . 0.000005 
ISSS^'^' 

U.OibS4_o.ooo34 

11.9 ±4.5 



3.5 -7.0'"' 
1 (54+0.19 

1.592851 ± 0.000045 
2026 ± 60 
n 02730+°"°'"'2 

U.UZ<OU„o. 00043 

3.0 ± 1.0 



4.8 ±0.8*''' 
1.68 ±0.09 
0.853585 ± 0.000024 
1800-2600 
0.0172 ± 0.00029 



1.5 



+0.8(d) 
-0.3 



(") Holman et al. <20T0t . 
('■^ Queloz et al. j2009t . 

f*^) Calculated value - assuming a Bond albedo of 0. 1 and a complete redistribution of heat for re-radiation, 
t'') , http://www.exoplanet.eu/, 



independently of the planet's composition, the hydrogen-helium gas envelope or water vapour atmosphere 
would escape within ^1 Ga, which is shorter than the calculated age of 1.2-2.3 Ga for CoRoT-7b. The age 
of Kepler-9 is 3±1 Ga and Kepler-10 is a relatively old star (11.9±4.5 Ga), so volatile-rich solutions are 
less likely if the Kepler-9 and Kepler-10 planetary systems are old and evaporation has been substantial. 
Jackson et al. ( 120101 ) studied the coupling of tidal evolution and evaporative mass loss on CoRoT-7b. Their 
investigation manifests that the orbital decay caused by the star-planet tidal interaction enhanced its mass 
loss rate. Such a large mass loss also suggests that, even if CoRoT-7b began as a gas giant planet, its original 
atmosphere has been completely evaporated now. As Kepler-9d and Kepler- 10b are all small size close-in 
planets (with orbital periods 1.59 and 0.84 days, respectively), evaporation on the planets is serious due 
to high stellar irradiation. Thus, the origin of low-mass exoplanets, like Kepler-9d and CoRoT-7b cannot 
be inferred from the present observations: they may have always had a rocky composition; they may be 
remnants of a Uranus-hke ice giant, or a gas giant with a small core that has been stripped of its gaseous 
envelope (Valencia et al. 120101 ). The latest results from Leitzinger et al. (1201 It indicated that hydrogen-rich 
gas giants within the mass domain of Jupiter and Saturn cannot thermally lose such an amount of mass that 
CoRoT-7b and Kepler-lOb would result in a rocky residue. They also concluded that these planets most 
Ukely were always rocky planets. 

In this paper, we assume that CoRoT-7b, Kepler-9d, and Kepler- 10b do not have significant gas en- 
velopes and consider them as fully differentiated with all possible iron-to-silicate ratios. The differentiation 
assumption is not restrictive since all the terrestrial planets and some large satellites in the solar system are 
known to be differentiated. Using a relatively simple but well suited model for terrestrial planets in solar 
system, we discuss the interior structures of CoRoT-7b, Kepler-9d and Kepler- 10b. We compare the silicate 
model and the carbon model (Seager et al. 12007) . Our main aims are: 1) infer their bulk composition and 
sharpen the observational constraints on their masses and radii, 2) compare their interior structures and infer 
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2 MODEL AND METHOD 

2.1 Model Equations and Numerical Method 

If the rotation rate is not extreme, self-gravitating bodies are close to spherical. Spherical structures are 
certainly an adequate starting place for studies of exoplanets. Three equations must be satisfied. They are 
the equation of hydrostatic equilibrium-Equation ([T), the mass conservation equation-Equation (|2]i and the 
equation of state (EOS)-Equation ^ for the material: 



dm (r) 



4Trr^p{r), (2) 



dr 

P(r)=/[p(r), r(r)]. (3) 

Where m (r) is the mass contained within radius r, P (r) is the pressure, T (r) is the temperature, and p (r) 
is the density of a spherical planet. The thermal contributions to the conclusion can be ignored which has 
been testified in Seager et al. (I2007l i and this approximation has been often used in other works (Fortney et 
al. 120071 Leger et al. l2004l l. As a reasonable approximation, the temperature-independent EOS is also used 
in this work. Using variable transition, if pressure is a nonlinear function of density. Equation ^ can be 
rewritten as, 

dp ^ m(r)p(r) I 

dr r2 ' dP{r)/dp{ry ^' 

We numerically integrate Equation (2) and Equation (4) starting at the planet's center r = and using the 
inner boundary condition m (0) =0 and p (0) = pc, where pc is a chosen central density. For the outer 
boundary condition, we use P (Rp) = 0. The choice of pc at the inner boundary and the outer boundary 
condition P (Rp) = define the planetary radius Rp and total mass Mp = m (Rp). Integrating Equation 
© and Equation (|4]i over and over for a range of pc provides the mass-radius relationship (see 3.1.1) for a 
given material. For differentiated model containing more than one kind of material, we specify the desired 
fractional mass of the core and of each shell. We then integrate Equation (|2]i and Equation ^ as specified 
above, given a pc and outer boundary condition. Using the planetary mass, we switch from one material 
to the next where the desired fractional mass is reached. Since we do not know the total mass that a given 
integration will yield ahead of time, we generally need to iterate a few times (change the value of pc) in 
order to produce a model with the desired distribution of material. 

2.2 The Choice of the EOS 

The materials used in this paper are ice VII, SiC, MgSi03 (enstatite, ab. en), (Mgo.ss. Feo.i2)Si03 (per- 
ovskite, ab. pv), Fe(a), Fe(e). Ice VII is used to model the ice composition in large moons of giant planets. 
MgSiOa and (Mgo.88> Feo.i2)Si03 are used to model the silicate mantles of terrestrial planets or large 
moons in our solar system. Enstatite is a low-pressure mineralogical phase of silicate, so it is only used for 
the validation of the model (see section 3.1.2). Fe(a) and Fe(e) are different phases of iron under different 
pressure, Fe(e) is the phase occurring at higher pressures. We compare their differences in mass-radius rela- 
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Table 2 Parameters for the Birch-Mumaghan EOS (BME) or Vinet EOS (VE) Fits. 



5 



Materials /\o (GPa) Kq po (kg/m^) EOS References 



ice VII 


23.7 ±0.9 


4.15 ±0.07 


1460 


BME 


1 


SiC 


227 ±3 


4.1 ±0.1 


3220 


BME 


2,3 


MgSiOa (en) 


125 


5 


3220 


BME 


2,4 


(Mgo.88, Feo.i2)Si03 (pv) 


266 ±6 


3.9 ±0.04 


4260 


BME 


2,5 


Fe(a) 


162.5 ± 5 


5.5 ±0.8 


7860 


BME 


2,6 


Fe(e) 


156.2 ± 1.8 


6.08 ±0.12 


8300 


VE 


7 



References: (1) Hemlev et al. ll987l (2) Ahrens l2000l (3) Alexsandrov et al. ll989l (4) OIinger ll977l (5) Knittle & 
JeanIoz [T987] (6) Takahashi & Spain [T989] (7) Anderson et al. [200T] 



cores. Let us consider the pressure scope in our topic. Pc, the central pressure of the planet can be estimated 
as following 

Pe ^ GA/2/4i?4. (5) 

It is about 10"^ GPa for the three exoplanets. For P < 200 GPa we use fits to experimental data, either 
the Vinet EOS (VE) Equation ^ or the Birch-Murnaghan EOS (BME) Equation In general, when 
P > 10^ GPa, where electron degeneracy pressure becomes increasingly important, the Thomas-Fermi- 
Dirac (TED) theoretical EOS (Salpeter & Zapolsky II967) must be employed, this is out of the range we 
need. The pressure range from approximately 200 to 10'^ GPa is not easily accessible to experiment nor is 
it well described by the TFD EOS (Seager et al. 120071 1. For all materials in this pressure regime we simply 
use the VE or BME. 

For a derivation of these EOSs, see Poirier (12000) . The VE is 



P = 3A-0 



2/3 



-1/3- 



cxp ■ 



(K'o 1) 



and the third-order finite strain BME is 



7/3 



5/3' 



1 ± ^ (A-'o - 4) 



2/3 



-1/3- 



(6) 



(7) 



Where po is the reference density, Kq — —V{dP/dV) is the bulk modulus of the material, Kq = 
(dK/dP) is its pressure derivative. Their values for the materials used in this paper are summarized in 
Table 2. The VE is used for the e-phase of iron because it is more suitable than the BME for extrapolation 
to high pressures (Poirier |2000t . The third-order BME is used for other materials. 

2.3 Carbon Model 



There is an assumption used in modeling the interior structure of solid exoplanets: their constituents are 
similar to the terrestrial planets in our solar system. They have an iron core and a silicate envelope for 
the differentiated model. But based on current observation, other possibilities can not be ruled out. The 
facts tell us that there are large differences between extrasolar planets and planets in our solar system. For 
example, research about hot Jupiter tells us it is danger to assume some planets can only form in particular 
area. It is useful to take as broad a view as possible (Fortney et al. 12007b . To other possible composition, 
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system were born in an environment where the carbon-to-oxygen ratio C/0 = 0.5. Their bulk compositions 
were determined by the high-temperature chemical equilibrium in protoplanetary disk. Silicates (i.e., Si-O 
compounds) are the dominant constituents in terrestrial mantles. But in an environment where C/O > 1, the 
condensation sequence changes dramatically (Larimer II 9751 Lewis [19741 Wood et al. 119931) . As a result, 
the high-temperature condensates will be carbon-rich compounds. Kuchner &. Seager (I2006I I have explored 
some planet-formation scenarios and they proposed that carbon planets composed largely of SiC or other 
carbides should form in such environment. Carbon-rich nebula created by the disruptions of carbon-rich 
stars or white dwarfs is a natural cradle for carbon planets, so planets around pulsars and white dwarfs are 
candidates of carbon planets. For example, the planets around pulsar PSR1257H- 12 (Wolszczan et al. 1 19921 1 
may be carbon-rich exoplanets. In the well-known /3 Pictoris debris disk (Roberge et al. l2006) . an exoplanet 
named (3 Pic b was discovered by Lagrange et al. ( I2008I I. it is most likely a carbon-rich planet also because 
it is located in a carbon-rich environment. Carbon planets can also form in a local area enriched in C or 
depleted in H2O in an otherwise solar abundance protoplanetary disk (Seager et al. 12007b . Lodders (I2004I I 
suggested that the planetary embryo that grew into Jupiter may have formed in such a locally carbon-rich 
area, and that Jupiter's embryo was a carbon planet. In this paper, we also consider the carbon model for 
CoRoT-7b, Kepler-9d, and Kepler- 10b, in which SiC is the major mantle constituent. The detailed formation 
scenarios for Earth-mass carbon planets can be learned from Kuchner &. Seager ( I2006I I. the latest research 
about the carbon-rich giant planets can be found in Madhusudhan et al. (1201 11 1. 

3 NUMERICAL RESULTS 

3.1 Validation in the Solar System 

3.1.1 Planets in the Mass-radius Relationship Diagram 

A mass-radius diagram is a useful tool when inferring a planet's bulk composition. Building on the work 
done by Zapolsky & Salpeter (11969b . we first consider planets of uniform composition. The positions of 
CoRoT-7b, Kepler-9d, and Kepler- 10b in mass-radius diagram can quickly tell us whether we can model 
them as solid planets. In Figure 1, the lines are mass-radius relationship curves for homogeneous planets. 
From top to bottom the homogeneous planets consist of ice VII, SiC, MgSiOa (enstatite, ab. en), (Mgo.ss; 
Feo.i2)Si03 (perovskite, ab. pv), Fe(a), Fe(£), respectively. Here, MgSiOa and (Mgo.ss. Feo.i2)Si03 are 
usually looked as representative material in mantle. The phase diagram of iron in Valencia et al. (2010) can 
help us to understand Fe(a) and Fe(£) phase of iron (under different pressure scope). CoRoT-7b, Kepler- 
9d, Kepler-lOb and the Earth are marked on the mass-radius diagram. Two notable exoplanets GJ 1214b 
(Charbonneau et al. 12009b and Kepler- lib (Lissauer et al. 1201 lb are introduced here as a comparison. Their 
measured masses and radii tell us that Kepler-9d, Kepler- 10b and CoRoT-7b are candidates of solid super- 
Earths. GJ 1214b is above the water ice VII line, this indicates it must contain some light elements such as 
H and He, most hkely it has a substantial gas envelope. Kepler- lib is above the solid SiC line, it should 
have a substantial water vapour or H/He gas envelope because it is so close to host star (10.3 days, Lissauer 
et al. l201 lb . As the validation, we also mark some planets and satellites (M < 1 M^) in our solar system 
in Figure 2. The Earth and the Venus fall between silicate and iron line, but close to silicate line, it implies 
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Table 3 Relative Error Between Calculated Radii and Actual Radii. 



Celestial Body Radius (km)'"' 


Assumed Composition'''' 


Model Radius (km) 


Relative Error (%) 


Mercury 


60%Fe + 40%en 


2474.43 


1.40 


2439.70 


60%Fe + 40%pv 


2344.78 


3.80 


Venus 


27.9%Fe + 72.1%en 


6283.54 


3.83 


6051.80 


27.9%Fe + 72.1%pv 


5862.68 


3.13 


Moon 


3%Fe + 97%en 


1745.05 


0.43 


1737.50 


3%Fe + 97%pv 


1652.37 


4.90 


Earth 


32.5%Fe + 67.5%en 


6589.91 


3.40 


6371.00 


32.5%Fe + 67.5%pv 


6269.81 


1.59 


Mars 


20%Fe + 80%en 


3409.08 


0.58 


3389.50 


20%Fe + 80%pv 


3149.51 


7.00 


Ganymede 


6.5%Fe + 48.5%en + 45%ice VII 


2525.19 


4.02 


2631.20 


6.5%Fe + 48.5%pv + 45%ice VII 


2454.14 


6.73 



|http://ssd.jpl.nasa.gov7| 



References: Mercury, Solomon <2003> : Venus, Zharkov J1983> : Earth, Ganymede, Seager et al. <2007> : Mars, 
Rivoldini et al. <20Tll . 



line, which indicates it must have significant amounts of ice composition. The Moon just lies on the silicate 
line, it agrees with the fact that the Moon is composed almost entirely of rock with a small iron core (Hu & 
Xu l2008l) . Unlike the Moon, Mercury must have a massive iron core. 

3.1.2 Code Test for the Differentiated Planets 

Mass-radius relationship diagram (or one layer model) can only give us a qualitative estimation, but it can 
not tell us the details of interior structure. To estimate the precision of our code, we calculate the radii 
of planets in solar system from their total mass and mass fraction published in relevant literatures. On 
observation, the planet mass and radius uncertainties are typically 5% to 10% (Selsis et al. 120071) . so a 
robust numerical code must guarantee the error is less than 5%. Table 3 are the results we get. We must 
emphasize that mass fractions of materials differ in literature, which indicates that there are still some 
uncertainties in the bulk composition of solar system objects. For example, the iron core mass fraction of 
the Moon is estimated at less than 10% of its total mass - 3% is found by Canup & Asphaug ( 1200 It . We 
use both pv and en phase of silicate to calculate the radii, the aim is to test the model approach only. Beside 
for Ganymede and Mars (both with the pv mantle), the relative error of all other planets and the Moon is 
less than 5%. Large relative error to Mars is understandable because Mars is relatively small, so pv is less 
suitable. Mercury is small too, but it has a large iron core, so pv or en can not cause large discrepancies. The 
relatively large error (6.73%) of Ganymede (with pv mantle) may be caused by the relatively low pressures 
in the deep interior of Ganymede (Sohl et al. I2002I I. and therefore pv is not very suitable, either The 
calculated Earth radius are within 101km (only 1.59% error) of the actual Earth, this means that pv is the 
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3.2 The Results for Three Exoplanets 

We examine the interior composition of Kepler- 10b, Kepler-9d and CoRoT-7b under the assumption of 
an iron core covered by a mantle composed of (Mgo.ss; Feo.i2)Si03 (silicate model, similar to Earth's 
mantle) or SiC (carbon model). When considering two-layer models, the measured mass and radius uniquely 
determine the planet's composition. The core mass fraction as a function of planet radius for the three 
planets are displayed in Figure 3-5. The solid red line denotes the fraction of planet's mass in its iron core 
according to an errorless planetary mass, while the green, white, and red shaded regions in Figure 3-4 
delimit the Icr, 2<t, and 3(t error bars on Alp, respectively. In Figure 5, the solid red line denotes the iron 
core mass fraction for ]\Ip = 5.25 (middle mass between 3.5-7.0 M^). The gray regions in Figure 5 
delimit the mass range we considered. The measured planet radius and its 1 cr error bars are denoted by solid 
blue and the dashed black vertical lines, respectively. 

Interior structure models can strengthen the observational constraints on a planet's mass and radius 
(Rogers et al. 12010b . For example, with the assumption that planets do not have a significant water or gas 
layer, some of the mass-radius pairs within Alp ± 1(Jm and Rp ± la^ (including the Oct mass-radius pair) 
can be ruled out because they correspond to bulk densities lower than a pure silicate planet. These excluded 
mass-radius pairs would necessitate water (or some other component lighter than silicate). In Figure 3-5, all 
mass-radius pairs in the la error bar (both mass and radius) compose an area. Here, we call them 'effective 
area'. 

To CoRoT-7b, if it has an Earth-like silicate mantle, its radius must be less than the measured value (1 .68 
i?0) and its mass has to be larger than the measured mass (4.8 M^) within Ict error bars, which implies an 
iron core mass fraction of less than ^ 27%. Within its observational value, its radius is unlikely less than 
the reliable value if it formed in the carbon rich environment and has a Moon-like composition (where the 
iron content is below 10% in carbon model). Its iron core is not heavier than ~ 70% of its total mass even 
at 3ct error bar (see carbon model in Figure 3), it is the largest value for both models. The silicate model 
of CoRoT-7b with an Earth-like iron core mass fraction (32.5%) is not consistent with the measured mass 
and radius within la. If its mass is ascertained with in 2a or 3ct error bar (in the range got by Hatzes et al. 
1201 II ), CoRoT-7b can be consistent with an Earth-like composition. From Figure 3, we can clearly see that 
the 'effective area' in the carbon model is much larger than for the silicate model. 

Kepler-lOb tells us other information. For its carbon model, there is a lower iron mass fraction limit 
(about 41%) within Ict error bar (see carbon model in Figure 4). KeplerlO (7.4-16.4 Ga) is older than 
CoRoT-7 (1.2-2.3 Ga), if CoRoT-7b is differentiated as we assumed, Kepler-lOb is more likely to be dif- 
ferentiated. This means that Kepler-lOb may have a larger iron core. Our simulation indicates Kepler-lOb 
is denser than CoRoT-7b, but the detailed calculation can tell us more: if CoRoT-7b has a silicate man- 
tle, whereas Kepler-lOb formed under carbon-rich circumstances, the upper core mass fraction limit of 
CoRoT-7b (~ 27%) is below the lower core mass fraction limit (41% mentioned above) of Kepler-lOb! 
To Kepler-lOb's carbon model, the upper limit of the iron core mass fraction is ^ 84% (see carbon model 
in Figure 4) within Ict error bar, this value is the largest in its two models and is larger as the suggested 
value of Mercury (^ 60%) - the densest planet in solar system. In Figure 6, we give the histograms of iron 



The silicate model and carbon rich model of CoRoT-7b, Kepler-9d and Kepler- 10b 9 

which we will discuss later. To Kepler- 10b, silicate model is compatible with an Earth-like composition 
within the Icr error bar, the mass is about 3.27 Af^ (lower limit of the observational value). The carbon 
model of Kepler-lOb is not compatible with an Earth-like core/mantle mass ratio (32.5% in core). 

Kepler-9d is shown in Figure 5, the gray band denotes the planetary mass range we are considering. If 
the maximum mantle-stripping theory is true, silicate model in Figure 5 tells us that Kepler-9d's radius can 
not be larger than ^1.75i?0. In detail, if it has an Earth-like core/mantle mass ratio, its radius is no more 
than ^-^1.7 Rg, with silicate model. Despite the mass of a transiting-only super-Earth is unknown, model 
simulations give constraints on their possible masses. For example, if Kepler-9d has an Earth-like mantle 
composition (silicate model in Figure 5), it has a lower mass limit according to its smallest possible radius, 
which is larger than the lower mass limit (3.5 Af®) discussed in Holman et al. (1201 Oi l. An Earth-like core 
mass fraction and the lower limit of observational radius (^1.50 i?®) yield a total mass of ^4.1 M^. In 
the mass range discussed in this paper, the carbon model of Kepler-9d is compatible with an Earth-like 
core/mantle mass ratio in a large scope of radius. 

Other information can be read from the Figure 3-5. For a given mass, the carbon model yields a larger 
radius than the silicate model. If the core is not pure iron but also contains a light element such as sulfur (in 
the Earth's core), all of the slantwise lines will sway to the right a little more (having no effect on the radius 
at a core mass fraction of 0.0). This indicates that the core mass fraction at a specified planetary radius 
will be larger Finally, we calculated the possible interior structure of CoRoT-7b and Kepler- 10b using their 
presumable masses (Figure 7-8). To CoRoT-7b, we take an Earth-like core/mantle mass ratio (core mass 
is 32.5%). From Figure 7, silicate model of CoRoT-7b can be ruled out because the calculated radius is 
outside the observational Icr error bar To Kepler-lOb, we take a Mercury-like core/mantle mass ratio (core 
mass is 60 %). In Figure 8, Kepler-lOb is compatible with a Mercury-Uke core/mantle mass ratio with its 
rehable mass for silicate model. 

4 DISCUSSION AND CONCLUSION 

Two recently discovered super-Earths, Kepler-9d and Kepler-lOb, and previously found CoRoT-7b are dis- 
cussed in this paper Because of their intense irradiation and small size, they are possibly absent of gas 
envelopes. Under the assumption of that they are mainly composed by refractory materials, their internal 
structures are deduced using suitable equation of state of the materials. The silicate and the carbon model 
are discussed in details. By analyzing the numerical results, we find some of the mass-radius pairs within 
Mp ± 1(Tm and Rp ± Icr/? can be ruled out. Therefore, our interior models sharpen the observational con- 
straints on their masses and radii. For CoRoT-7b, an Earth-like composition with a silicate (pv) mantle is 
less likely within both the mass and the radius Icr error bars, but a Mars-like composition may be suitable. 
To Kepler-lOb's two models, a Mercury-like composition having 60% of its mass in iron core is consis- 
tent with the measured mass and radius within 1 cr error bar There is an upper and lower limit for the iron 
core mass fraction of every planet. To CoRoT-7b, the upper iron core mass fraction limits are 27% (silicate 
model) and 56% (carbon model). For Kepler-lOb, the silicate model have an upper limit of about 75% and 
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9d can also be studied, for example, its radius is no more than 1.7 i?^ in silicate model. The lower mass 
limit (3.5 Af® we adopted) is also ruled out in its silicate model. 

Our quantitative calculation tells us Kepler- 10b may be a Mercury-like planet. CoRoT-7b and Kepler- 
10b share the similar features: their radii, their period, the type of host star and perhaps their mass. Why 
is Kepler-lOb denser than CoRoT-7b? Their different ages may give us some answers. To the anomalously 
dense Mercury in solar system, there are three competing viewpoints (Solomon l2003l l: 1) the solar nebula 
caused aerodynamic drag on the particles from which Mercury was accreting, the iron and silicate have 
different response to this drag, so the lighter particles were lost from the accreting material at the onset of 
accretion. 2) it surface rock could have been vaporized by strong radiation from a hot nebula and have been 
carried away by the solar wind, 3) much of the original crust and mantle have been stripped away by a giant 
impact caused by a planetesimal. Because of lacking detailed observational data, which hypothesis is true 
can not be determined. However, each hypothesis predicts a different surface composition. So two promising 
space missions, underway MESSENGER (Solomon et al. 120011 NASA) and the upcoming BepiColombo 
(Grard et al. 120011 ESA), could bring us the answer The second and third explanations invoke processes 
late in the planetary accretion process, after the protoplanet had differentiated silicate mantle from metal 
core. If CoRoT-7b and Kepler-lOb were bom in a similar environment, Kepler-lOb's abnormal iron content 
can be naturally associated with its older age. We think the postnatal alteration is more likely. 

The lower limit of iron content of a super-Earth is an interesting question because a pure silicate super- 
Earth (iron mass fraction is zero) is unimaginable under differentiation hypothesis. Under physically plau- 
sible conditions analogous to solar system, planet formation will lead to differentiated super-Earth with an 
iron core covered by a silicate mantle, with the proportions of each determined by the local Si/Fe ratio 
(Grasset et al. |2009l l. The only way to significantly increase the mean density of a planet requires removal 
of an extended part of its silicate mantle. Marcus et al. (2010) proposed an efficient method - giant impacts 
between super-Earths. On the other hand, they prefer a lower limit of iron content of 33% by considering 
the standard cosmic abundances (see also Valencia et al. 12007b . The initial Fe/Si ratio (average value) is 
used to give a lower limit for the iron content of super-Earth. But, in practice, formation progress must be 
influenced by other factors (such as the survival competition between planets) because some planets in solar 
system have an iron content lower than this 'lower limit' . We think that referring to the iron content in chon- 
drites is a sound choice. L chondrites have lower total iron content, which is about 10% (Chen et al. |19941 l. 
L chondrites are the largest group in the ordinary chondrites (40 % , pittp : //w w w. nhm. ac .uk^ . Therefore, we 
use 10% as the lower limit for silicate model. Regarding Figure 6, if there is a calculated lower Umit of the 
iron content, we will use it, otherwise we use 10% (only for silicate model). 

Our main conclusions are summarized as follows: 1) The iron core of CoRoT-7b is not more than 27% 
of its total mass within 1 a mass-radius error bars, so an Earth-like composition is less likely, but its carbon 
rich model can be compatible with an Earth-like core/mantle mass fraction; 2) Kepler-lOb is more likely 
with a Mercury-like composition, its old age implies that its high iron content may be a result of strong solar 
wind or giant impact; 3) the transiting-only super-Earth Kepler-9d is also discussed. Combining its possible 
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Kepler-9d. The radius derived from transit method and the mass got by radial velocity detections together 
with the model simulation can constrain the possible compositions of exoplanets. 

Credible conclusions desiderate comprehensive observations. Planetary formation theories, thermal evo- 
lution models, and studies of the cosmic abundance of elements can be used to place additional constraints 
on a planet's interior composition. Besides the temperature-independent model used in this paper, more 
complicated models of solid exoplanets can be found in other literatures. For example, Valencia et al. 
(12006b and Sotin et al. (12007b considered the temperature profile and mixture of materials in the mantle or 
core. Wagner et al. (1201 lb used material laws in the infinite pressure limit to improve the model of solid 
exoplanets. In the near future, high-sensitivity spectroscopic transit observations of these exoplanets should 
constrain the compositions of the evaporating flow and therefore allow us to distinguish between rocky or 
gas-rich planets. In particular for exoplanets, without any opportunity to in situ composition measurements 
and gravitational moment measurements from spacecraft flybys, we will be permanently limited in what we 
can infer about the interior composition from their observational mass and radius (Rogers &. Seager l2010b . 
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Fig. 1 Mass-radius relationship diagram for Mp < 10 M^. The solid or dashed lines are ho- 
mogeneous planets. Five exoplanets and the Earth are marked on the diagram. The mass range 
of Kepler-9d (red error bar) is taken from Holman et al. (1201 Oi l. The locations of CoRoT-7b, 
Kepler-9d and Kepler-lOb indicate they can be described as solid planets. 
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Fig. 2 Mass-radius relationship diagram for Mp < 1 Af^. The solid or dashed lines are ho- 
mogeneous planets. Some planets and satellites in solar system are shown in it. Their major 
constituents can be estimated using this mass-radius diagram. 
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Fig. 3 Core mass fraction as a function of planetary radius for CoRoT-7b. Silicate model 
((Mgo.88, Feo.i2)Si03) and carbon model (SiC) are discussed. The planetary mass is M = 
4.8 ± 0.8 Mq. The green, white and red regions denote the core mass fractions obtained by 
varying its mass within la, la and 3a error bars. The dashed black vertical lines delimit the 
measured radius R = 1.68 ± 0.09 i?©. The Silicate model of CoRoT-7b had been discussed by 
Rogers & Seager (1201 Oi l and Valencia et al. (1201 Oi l. We redo it for the comparison with its carbon 
model. 
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Fig. 4 Kepler-lOb core mass fraction as a function of planetary radius. Silicate model ((Mgo.ss. 
Feo.i2)Si03) and carbon model (SiC) are discussed. The lines and areas are similar to Fig. 3. 
The planetary mass is M — 4.56l|^'2g ^^ffi '^^e measured radius is i? = 1.416lQQ3g -^ffi- 
/S.M = ±1.29 Afg is used for la and 3<t error bars of its mass. 
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Fig. 5 Kepler-9d core mass fraction as a function of planetary radius. Silicate model ((Mgo.ss. 
Feo.i2)Si03) and carbon model (SiC) are discussed. The measured radius is i? = 1.64^0 '^[4 Rq. 
The gray region denotes the core mass fractions obtained by varying the mass within the mass 
range we considered, Af = 5.25 ± 1.75 Mq. 
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Fig. 6 The histograms of iron content for CoRoT-7b, Kepler- 10b and other celestial bodies men- 
tioned in our work. The corresponding references are the same in Table 3. The green part denotes 
the upper and lower limit for silicate model, the blue part denotes similar limit for carbon model 
(see details in Section 4). To the grid range on the histogram of Kepler-lOb (silicate model), the 
lower limit (33%) adopted by Marcus et al. ( 1201 Oi l is used. Obviously, from a statistical point of 
view, CoRoT-7b is compatible with a Mars-like core/mantle mass fraction, whereas Kepler- 10b 
is compatible with a Mercury-like composition. 
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Fig. 7 The interior structure of CoRoT-7b. We use the reliable mass AI ~ 4.8 with an Earth- 
like core/mantle mass ratio (iron core is 32.5% ). Two mantle compositions are considered, the 
red line is for (Mgo.ss. Feo.i2)Si03 and the blue line is for SiC. Its radius scope R = 1.68 ± 
0.09 i?0 is shown using vertical lines. Silicate model can be ruled out because the corresponding 
radius is out the 1 cr error bar for its reliable mass. 
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Fig. 8 The interior structure of Kepler- 10b. We use a Mercury-like core/mantle mass ratio (iron 
core is 60% ) and its reliable mass M = 4.56 A/®. Two mantle compositions are consid- 
ered, the red line is for (Mgg.ss. Feo.i2)Si03 and the blue line is for SiC. Its radius scope 
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is shown using vertical lines. Carbon model can be ruled out because 



the corresponding radius is out the 1 a error bar for its reliable mass. 



